Role of surface anisotropy for magnetic impurities in electron dephasing and energy 

relaxation and their size effect 
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Recently the electron dephasing and energy relaxation due to different magnetic impurities have 
been extensively investigated experimentally in thin wires and in this Letter these quantities are 
theoretically studied, ft was shown earlier that a magnetic impurity in a metallic host with strong 
spin-orbit interaction experiences a surface anisotropy of the form H = Kd(nS) 2 which causes size 
effects for impurities with integer spin. Here we show that the dephasing and the energy relaxation 
are influenced by the surface anisotropy in very different ways for integer spin having a singlet 
ground state. That must result also in strong size effects and may resolve the puzzle between the 
concentrations estimated from the two kind of experiments. 
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In the present Letter, to our knowledge, is the first time 
a mechanism is presented which can resolve the seriously 
puzzling observation that in some cases the influence of 
magnetic impurities on electron dephasing and energy 
relaxation have drastically different strengths. The sur- 
prising difference in the estimated impurity concentration 
raised the doubt about the role and even the presence of 
magnetic impurities [lj. That mechanism is based on 
surface magnetic anisotropy 2] which is resulting in a 
strong size dependence. 

The size dependence of the Kondo effect |j3, Ll was dis- 
covered more than ten years ago and since than it has 
been carefully studied experimentally @, IE H S lill • 
That cannot be attributed to the size of the Kondo 
screening cloud reduced by the size of the sample as 
only the energy separation of the metallic electron levels 
are relevant. That problem was resolved by the sugges- 
tion that the magnetic impurities in the metallic host 
with strong spin-orbit interaction experience a surface 
anisotropy |2(. 

In mesoscopic metallic systems the electron dephasing 
and energy relaxation are the central issues in under- 
standing their transport properties [llj . The interest has 
been intensified by the debate over the saturation of de- 
phasing at low temperature |l2j . The dephasing is deter- 
mined e.g. from measurements of magnetoresistance and 
Aharonov-Bohm rings in magnetic field |l2j, |ljj, I hi Il5| 
while the energy relaxation from transport in short wires 
is found by determining the nonequilibrium electron en- 
ergy distribution as shown by the Saclay group |lq|. In 
many cases the deviations from the expectations of the 
theory of Altshuler, Aronov and coworkers [ljj were at- 
tributed to the presence of magnetic impurities either 
implanted or contained by the starting material as con- 
tamination. In addition to the energy relaxation due to 
electron-electron, clcctron-phonon interaction the mag- 
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FIG. 1: The magnetic impurity at distance d from the sur- 
face. 



netic impurity mediated electron-electron interaction can 
play an essential role, which is supported also by new ex- 
periments in magnetic field |18| . It has been known since 
a long time |l9| that such interaction is singular in the 
energy transfer E and recently Kaminski and Glazman 
[2JJ| called the attention to similar 1/E 2 singularity in 
the electron-electron scattering rate phenomenologically 
suggested by the Saclay group. Using that mechanism 
the electron transport is determined by the Boltzmann 
equation and compared with the experimentally deter- 
mined electron distributions and the impurity concentra- 
tions were adjusted [J |2j, |22, |23|, |24( • In some cases the 
estimated magnetic impurity concentrations using that 
method are much larger than those determined from the 
dephasing rate, even by two orders of magnitude. In 
AuPd samples with large spin-orbit interaction the size 
dependence of dephasing was also observed [14j where 
for smaller size the electron-electron interaction domi- 
nates while for larger samples there are additional scat- 
tering mechanism resulting in saturation. That trend is 
just the opposite what could be expected in case of addi- 
tional scattering centers at the surface. In the following 
the possible role of surface anisotropy in these phenom- 
ena is discussed. 

Surface anisotropy: 



For plane-like surfaces (see Fig 
surface anisotropy has the form 



i;2] 



the spin-orbit-induced 



H = A' d (nS) 



(1) 



where n is the normal direction of the experienced sur- 
face element and S is the spin of the impurity. The 
anisotropy constant K d > is inversely proportional to 
the distance d measured from the surface thus it has 
the form K d = a 4. For thin films using the assump- 
tion that the to surfaces of the film contribute additively, 



K d 



t~d> 



a was obtained from fitting the Kondo 



resistivity of Au(Fe) and Cu(Fe) films as a = 247AK 
25f. of magnetoresistance of Au(Fe) films as a = 42AK 
26], and from multilayer experiments on Au(Fe) films as 
a = 60AK 7]. The parameter a depends also on the 
disorder on the surface and in the bulk. 

According to the anisotropy there are different split- 
ting schemes for integer and half- integer spins. For inte- 
ger spins (e.g. Fe, Cr S = 2) the ground state is a singlet, 
whereas for half-integer spins (e.g. Mn S = 5/2) it is a 
Kramers doublet. Thus for integer spins the anisotropy 
causes size effects e.g. in Kondo resistivity |3J, magne- 
toresistance [3, thermopower J6j, impurity spin magne- 
tization 8] , but for half-integer spins not [J, |jj, |27J • It is 
demonstrated that there is a crucial difference between 
the cases of integer and half-integer spin. That difference 
can be less pronounced for e.g. S — 5/2 as in the spin 
glass region pairs or clusters can be formed which could 
have also integer spins showing size dependence [lOj. 
Dephasing: 

As the experiments are carried out at low temperature 
e.g. T ~ 40mK, thus in thermal equilibrium most of 
the higher levels cannot contribute to any dynamics for 
impurities with large enough anisotropy, K d > kT. In 
case of integer spins the impurity is frozen in a singlet 
ground state which cannot lead to dephasing, in contrary 
to the half-integer spin case where the lowest states form 
a doublet and we do see dephasing. Samples with Fe 
implantation or contamination must show a strong size 
dependence in contrary to Mn. 
Energy relaxation: 

The nonequilibrium distribution function of a metallic 
wire with length L and bias U in the diffusive limit is 
determined by the Boltzmann equation 



df(e,x) 
dt 
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fdE{f(e)[l-f(s-E)]W(e,E) 
[l-f(e)]f(e-E)W(s-E,-E)} (2) 



where W(e,E) is the scattering rate, td = 75- is the 
diffusion constant, / is assumed not depending on the 
spin, and x denotes the position in the wire in the units of 
L. Starting with the solution without inelastic scattering 
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FIG. 2: The diagrams used for calculating (a) the kernel and 
(b) the Korringa lifetime of the impurity spin. The solid lines 
denote the conduction electrons, the dotted lines the impurity 
spin, and the blob is the Kondo coupling. 



mechanism 
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/ w (£,z) = (1 - x)n F (e - — ) + xn F (e + — ) (3) 

and taking into account inelastic scattering in W , the 
Boltzmann equation can be solved self-consistently at 
least numerically. 

Here we examine the effect of the surface anisotropy on 
the energy relaxation. Similar to the case of finite mag- 
netic field [l| the first order processes contribute also to 
the scattering rate and the spin occupation numbers PmS 
depend also on the voltage U. Calculating them from the 
first order processes we solved the Boltzmann equation 
self-consistently using the following collision integral 

Icl ({/}) =/ dE J de'K s MM , (E, e, e', K d ){p M f(e)f(s') 

x[l-f{e-E)][l-f(e' + E + K d M 2 - K d M' 2 )} 

-p M ,[l-f(e)][l-f{e')]f{s-E) 

xf(e' + E + K d M 2 - K d M' 2 )} (4) 

where the kernel Kfj M , describes electron-electron in- 
teraction mediated by Kondo impurities with surface 
anisotropy. 

For simplicity we considered the S = 1 case when the 
x-dependent pms are determined from the first order pro- 
cesses by 



Po _ J defje, x)(l - /(e + K d , x)) 
pi fdef(s,x)(l-f(E-K d ,x)) 



(5) 



and 2pi + po = 1. 

The kernel in Eq. (0J can be calculated from the di- 
agram Fig. |2 (a) in the Kondo model with anisotropy 

H = Y, £ ka{ a a ka + K d {nS) 2 

k,a 

+ 2^ Jmm'Smm' (a k(j cr aa >ak'a'), (6) 
MM' 



where a\ a (afe CT ) creates (annihilates) a conduction elec- 
tron with momentum k, spin a and energy Ek measured 
from the Fermi level, a stands for the Pauli matrices and 
Jmm''s are the Kondo couplings. The dependence of the 



interaction kernel Kf [M , on the energy transfer E for 
t k = oo is {E 2 )-\ (E±K d )~ 2 , (E + Ka)- 1 ■ (E~K d )- 1 
in different terms, respectively. 

For sake of simplicity we used an appropriate constant 
value J instead of the renormalized Kondo couplings de- 
pending on M, M' . The influence of such an approxima- 
tion was examined in a preceding self-consistent calcula- 
tion without surface anisotropy |28j . There the renormal- 
ized coupling was calculated as the solution of the leading 
logarithmic scaling equation assuming similar resumma- 
tion as in equilibrium and smeared by the spin spectral 

function with finite Korringa lifetime p s (s) — i - 2TK h i . 

The validity of the logarithmic approximation was always 
checked by plotting the actual Kondo coupling. From the 
numerical calculations we can conclude [28| that to get 
the same results it is a good approximation to replace the 
renormalized coupling in the kernel by an appropriately 
chosen constant value. Furthermore, the smearing of the 
renormalized coupling has very small effect |28| on the 
results for the parameters consistent with the experimen- 
tal situation and our results were in complete agreement 
with Ref.Q]. 

As the weak dependence on the Korringa lifetime |29j 
tk of the impurity spin we used the value for K d = 
calculating it from the diagram Fig. (b) as 

2tt( Po J) 2 S(S+1) fde(l-f(e,x))f(e,x). 



strength of the anisotropy constant K d is illustrated in 
Fig. |21 Increasing K d first the energy transfer is get- 
ting larger but for larger K d the ground state is frozen 
in, similar to the magnetic field dependence discussed in 
Ref. 0] . We can conclude that the contribution of mag- 
netic impurities is enhanced or unchanged in case of finite 
anisotropy K d < ell. For K d ~ 0.1— 0.2K which is a good 
estimation for the strength of the anisotropy for the wires 
with width of ~ 45nm and thickness of ~ 85— HOnm used 
in the experiments, the energy relaxation is only slightly 
affected by the anisotropy. 




2t k {x / 

(7) 
where po is the conduction electron density of states for 
one spin direction. 

At each step of the iteration solving the Boltzmann 
equation self-consistently, both the spin occupation num- 
bers and the Korringa lifetime were calculated from the 
actual /. 



FIG. 4: Fit on the experimental data of Cu wires at x = 0.485 
[l8| by the calculated distribution function for different Kd 
and c pairs. The other parameters are U = O.lmV, poJ — 0.11 
and td = 2.8ns. 

The goal of the present Letter is not to make opti- 
mal fitting of the experimental curves and determine the 
value of K d which must have a broad distribution itself. 
We demonstrate, however, that the experimental curves 
can be fitted by using combinations of different values of 
the concentration c and of K dl and the larger the c is the 
smaller the necessary K d . As a demonstration we com- 
pare our results to the experimental data on Cu wires 
at x = 0.485 [1S| in Fig. 0] as for Cu wires the impuri- 
ties may be CuO on the surface having S = 1 spin [30j. 
The other fixed parameters are U = O.lmV, poJ = 0.11 
and td — 2.8ns, and similarly good fits are obtained for 
U = 0.3mV as well. The parameters are somewhat dif- 
ferent, which is not surprising as the distribution for K d 
is not taken into account. It is important to note, that 
in some cases the origin of the magnetic impurities is not 
known, therefore the Kondo temperature corresponding 
to poJ in our simple approximation is also a fit parame- 
ter. 

The half-integer case must be very similar to the case 
without surface anisotropy because of the degeneracy, 
and only the spin dependent prefactors are different. 

The two-level system (TLS) may result in somewhat 
similar behavior |31j. If T << A, where A is the split- 
Trie dependence of the distribution function on the ting, the dephasing is blocked. In the nonequilibrium 




FIG. 3: The calculated distribution function at x — 0.485 for 
different strength of the anisotropy constant Kd- The other 
parameters are U = O.lmV, c = 8ppm, poJ — 0.11, and 

td = 2.8ns. 



case with applied voltage U (A < ell ) th e spin dynamics 
reenters and could lead to dephasing 32] similarly to the 
anisotropy case. Similarly, the energy relaxation becomes 
also possible but to get -^ singularity at least two non- 
commuting couplings are needed [33|, thus interaction 
describing electron screening and electron induced tran- 
sition between the levels are required |3J| . In this case the 
splitting must be small A < ell, but the coupling can be 
weak enough to be outside the Kondo region. That may 
result in weak, magnetic field independent contribution, 
what is suggested by the experiments |35| . 

In summary, the surface anisotropy for integer spins is 
suggested to reduce drastically the dephasing rate, while 
the energy relaxation is much less influenced. In the first 
case for low temperature and thermal equilibrium the 
spin dynamics and therefore the dephasing are frozen out 
while in the out-of equilibrium metallic wire experiments 
that can reenter. That suggests a pronounced size depen- 
dence and very different concentration for the dynami- 
cally active impurities in the dephasing and the out-of- 
cquilibrium wire experiments. In the case of half-integer 
spin having a Kramers doublet as the lowest state instead 
of a singlet these cannot be expected. Further careful 
experiments for the size dependence and implanted im- 
purities are required. 
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